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Abstract 

A rigorous reduction of the many-body wave scattering problem to 
solving a linear algebraic system is given bypassing solving the usual system 
of integral equation. The limiting case of infinitely many small particles 
embedded into a medium is considered and the limiting equation for the 
field in the medium is derived. The impedance boundary conditions are 
imposed on the boundaries of small bodies. The case of Neumann boundary 
conditions (acoustically hard particles) is also considered. Applications 
to creating materials with a desired refraction coefficient are given. It 
is proved that by embedding suitable number of small particles per unit 
volume of the original material with suitable boundary impedances one 
can create a new material with any desired refraction coefficient. The 
governing equation is a scalar Helmholtz equation, which one obtains by 
Fourier transforming the wave equation. 



1 Introduction 

This paper is a continuation of [H] and uses some of the results from [9], [7j, 
[TO] . [T5] . [T6] . Applications of our theory to creating materials with desired 
refraction coefficient, including negative refraction, are discussed in [11], |12j . 
[13], [17]. Wave scattering by small bodies is a classical branch of science: it was 
originated by Rayleigh in 1871. In [2] one finds a discussion of wave scattering by 
a small particle. In [I] there is a review of the low frequency scattering theory and 
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formulas for scattering by small balls and ellipsoids are given. In [9] the theory is 
developed for small bodies of arbitrary shapes. In [13] the many-body scattering 
problem was reduced to solving linear algebraic systems bypassing the usual study 
of a system of integral equations. In this paper we apply the approach proposed 
in [13] and study the limiting behavior of the scattering solution when the number 
of small bodies tends to infinity in such a way that the characteristic size a of the 
small particles is related to their number M so that M = O(-) in Theorem 2, and 
M = O(-) in Theorem 3. Sufficient conditions for convergence of the scattering 
solution in this limiting process are given. We prove that these conditions are, in 
some sense, also necessary for convergence. The limit of the scattering solution 
is a function, which satisfies some differential or integral-differential equations. 
These equations describe the behavior of the wave field in the new medium, 
obtained in the limit. 

There is a large literature on the calculation of the effective dielectric permit- 
tivity and magnetic permeability of the composite materials (Maxwell-Garnett 
and Bruggeman recipes and their numerous versions, and newer theories [18J, 
[3]). In the literature mostly a randomly uniform distribution of the inclusions 
is assumed and the resulting homogenized medium is described by effective con- 
stant dielectric permittivity and magnetic permeability, which can be tensors. In 
this work the propagation and scattering of scalar waves are discussed, and the 
"homogenized" medium is described not by a constant refraction coefficient, but 
by a refraction coefficient which is a function of spatial variables. 

Let us formulate the problem. Consider first a bounded domain D C M 3 
filled with a material with a known refraction coefficient n (x). The governing 
equation is: 

L u := (V 2 + k 2 n (x))u = in M 3 . (1.1) 

We assume that n (x) = 1 in D' — M?\D, k = const > 0, and no = 
m&x x£D \n (x)\ < oo. The operator L can be written as a Schrodinger oper- 
ator: 

L = \7 2 + k 2 -q (x), q (x):=k 2 [l-n (x)}, (1.2) 
and go = in D' . One has 

no (a:) = 1 - k~ 2 q (x), 

so there is a one-to-one correspondence between n (x) and qo(x). If n (x) is 
known, then one knows the scattering solution: 

L u = in R 3 , 

u (x) =e ika - x + A ((3 } a) — + o(-), r = \x\ -»• oo, /3 := -. (1.3) 

The coefficient Aq(P, a) is called the scattering amplitude, the unit vector a G S 2 
is given, a is the direction of the incident plane wave e lka ' x , S 2 is the unit sphere 
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in M 3 , (3 G S 2 is the direction of the scattered wave, k > is a wave number, 
which we assume fixed throughout the paper. By this reason we do not show the 
^-dependence of A and uq. 

Let G(x, y) be the resolvent kernel of Lq satisfying the radiation condition (or 
the limiting absorption principle): 

L G(x,y) = -5(x -y) in M 3 . (1.4) 

This function G(x,y) is known because qo(%) is known. 

Consider now the scattering problem for many small bodies D m embedded in 
D, 1 < m < M: 

M 

L u M = inM 3 \ |J D m , (1.5) 

m=l 

um = Uo + A M (f3,a) ho(-), r = \x\ — >• oo, - = (3, (1.6) 

-Qj^r = CrnU M On S m '.= dD m , 1 < 771 < M, (1.7) 

where Uq is the solution ot the scattering problem (ll.3p . Here N is the normal 
to S m pointing out of D m , ( m is a complex number, the boundary impedace, 
Im(m < 0, S m is uniformly C 1,A with respect to m, 1 < m < M. By C 1,A surface 
we mean the surface with local equation x 3 = f(xi,X2), where / G C 1,x , A > 0. 
We assume throughout this paper that 

noka <C 1, d^> a, (1.8) 

where ^ 

a = - max diamD m , <i = min dist(D m , DA (1.9) 

By Kn := \D m \ the volume of D m is denoted, and by \S m \ the surface area of S m 
is denoted. 

One can prove (see Section [3]) that problem (11. 5p - (11. 7\i has at most one 
solution if lm( m < 0, 1 < m < M, and Img (^) < 0. 

We look for the solution to problem (11. 5p - (ll.7p of the form 

M 

u M {x) = u (x) + ^ / G{x,s)a m {s)ds, (1.10) 

m=l 

where a m should be found from the boundary conditions (\1.7\i . For any a m the 
function fll.lOp solves equation (11.51) and satisfies condition (11.61) : 

i M r 

A M ((3,a) = —22 u (s,-p)a m (s)ds. (1.11) 

m=l J Sm 
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Formula (11. lip follows from (11.61) . (11.101) and the Ramm's lemma ([8], formulas 
(5.1.31), (5.1.36)): 

gifc|a;| / 1 \ X 

G(x, y) = - — : — - u (y, a) + o( - — - J , \x\ ->• oo, -p- r = -a, (1.12) 
47r|x| V|x|/ \x\ 

where Uq(x, a) is the scattering solution. A similar formula was proved earlier in 
[Sj, P- 46, for the resolvent kernel of the Laplacian in the exterior of a bounded 
obstacle, (and even earlier, in [5], for some unbounded obstacles). The scattering 
amplitude for problem (11.51) - (11.71) is 

A(J3,a) = Ao(p,a)+A M (J3,a), (1.13) 

where A is defined in ( 11. 3ft and A m is defined in ( II. 61) . If ka is sufficiently 
small, then k 2 is not a Dirichlet eigenvalue of the operator V 2 — qo{x) in D m , 
1 < m < M. If 

ImCm<0, l<m<M; Img (x) < 0, (1.14) 
then the unique solution to problem (11. 5p - (I1.7P can be found in the form (ll.lOp . 

Theorem 1 Assume (11.81) and (11.141) . Then problem (11. 5p - (11.71) has a solution 
of the form (ll.lOp and this is the unique solution of the problem (11. 5p - (11.71) . 

Proof of Theorem [1] is given in Section [3j In [3] there is a detailed study of 
boundary value problems in domains of the type D' := R 3 \(J^ =1 D m . In [3] 
the case of Dirichlet boundary condition on S m was studied, the case of Neu- 
mann boundary condition was mentioned as an open problem, and the case of 
impedance boundary condition was not studied. 
Let 

e ik\x— y\ J 

9(x,y) ■=-—, r, g {x,y) := — r. (1.15) 

Air\x — y\ Air\x — y\ 

Note that 

G(x,y) = g(x,y) - [ g(x, z)q (z)G(z,y)dz. (1.16) 
We need two lemmas. 
Lemma 1 // 



D 



then 



t — x\ < a, \x — y\ > d 3> a, (1-17) 

ka s 



\g(t,y)-g(x,y)\<c(^ + -^), (1.18) 
where c > stands for various positive constants independent of a and d. 
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Lemma 2 If (11.171) holds, then 

\G(t,y)-G(x,y)\<c^ + ^-). (1.19) 

These lemmas are proved in Section [3j 

Let us formulate our results under simplifying but physically reasonable as- 
sumptions. 

Theorem 2 Assume that 

lim ^ m m = h(x m ), where J m := [ [ -, (1.20) 



a^o 47r|5' m | ' "' lq Iq \s — t 



and /or am/ subdomain D G D the following relation holds 

y i = I / Ar(x)dx. (i.2i) 

Assume that \S m \ = cia 2 , and J m = C20? , where C\, C2 > are constants indepen- 
dent of m. 

Finally we assume that M = 0(a~ x ) and d = 0(a^ 3 ) as a — > 0. Under these 
assumptions there exists the limit: 

lim um(x) = u(x) = u(x,a). (1-22) 

M^oo 

This u(x) solves the equations: 

u(x)=u Q (x)- [ G(x,y)p(y)u(y)dy, 

and 



D 



Lu := [V 2 + k 2 - q(x)}u = m R 3 , (1.23) 
where the potential q is of the form: 

1 \ / \ / \ 1 \ 47tc? N(x)h(x) 

g(x) = q Q (x) + p i , p i = * , 77 / , (1-24) 

c 2 [l + n(x)J 

and u satisfies the radiation condition: 

ikr . 1 . 

u = e ika ' x + A(P, a) + o(-J, r = b| -> 00, (1.25) 

r Vr / 

w/iere 

A(f3,a) = A (f3,a) - ^ J u (y, -ff)p{y)u{y,a)dy, (1.26) 
and 1*0(2/, ~~ /?) scattering solution defined in ( II. 3p . 



Theorem 3 Assume that ( m = 0, 1 < m < M , and the following limits exist: 

lim n E = I Ppi(v>(v)*v> i 1 - 27 ) 

D m cD 

lim V V m = f u(y)dy, (1.28) 



D m CD 



where u(y) > and (3 P j(y) are continuous functions in D, and (3^ is the mag- 
netic polarizability tensor of the body D m , defined in (12.381) -( 12. 391) . see below. 
Then the function um{x), defined in (11.101) . tends to the limit: 

lim um(x) = U{x) = U(x, a), (1-29) 



andlA{x) solves the equation: 



+ / G(x,y)AU(y)v(y)dy- V / 
Jd P,i=i Jc 



dG(x,y) dU{y) 



U(x)=u (x) + / G(x,y)AU(y)u(y)dy - \ / — —2 [3 pj (y)v(y)dy. 



: , 3 zJd dy p d Vj rm 

(1.30) 

If all the small particles are balls of radius a > 0, then 

V m = -^r - , \S m \ = 47ra 2 , J m = 167r 2 a 3 , / ^ = a, s G £> m . 
3 /o 47T s — t 



In this case 



dt Iff dtds 



s I $ t\ Srn \ I <> I Q 47T \S t 



that is, the mean value of the integral J s r^r on the surface S m equals to this 
integral. If S m is not a sphere, this mean value is an approximate value of the 
above integral. 

Note that under the assumptions of Theorem [2] one has M = 0(a^ 1 ), while 
under the assumptions of Theorem [3] one has M = 0(a~ 3 ) (see formula (I2.50P be- 
low). Therefore, one needs many more particles to deal with the Neumann bound- 
ary condition, that is, with acoustically hard particles, than with the impedance 
boundary condition with large boundary impedance ( = 0(a _1 ). We will discuss 
at the end of Section H] in more detail the question concerning the compatibility 
of the assumption (11. 8B . namely d ^> a, and the existence of the limits ( 11.271) 
and (11.281) . It will be shown that the assumption d ^> a is compatible with the 
existence of the limit (11.281) only if v{y) is sufficiently small, and in this case the 
existence of the limit (I1.27P is also compatible with the assumption d^> a. 

In Section [2] Theorems [2] and [3] are proved. In Section [3] Theorem [3] and 
Lemmas 1, 2 are proved. In Section H] some examples are given, the significance 
of the compatibility of the assumptions d ^> a and (11.211) . (11.271) - (11.281) is 
discussed, and a possible application of our results to creating materials with a 
desired refraction coefficient is described. 
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2 Proof of Theorem [2 



Let us look for the solution to problem ( 11. 5p - ( 11.7j) of the form: 

M 

um = u {x) + ^2 G(x,s)a m (s)ds, (2.1) 

m=l ** Sm 



where G(x,y) is the resolvent kernel of L , see (jl.4p . and arbitrary func- 

tions at the moment. For any cr m the function (12.11) solves equation (II. 5p and 
satisfies the radiation condition (jl.fip . Since problem (11.51) - (jl.7p has at most 
one solution, the function (12.11) is the unique solution to (11.51) - (II. 7p provided 
that a m are chosen so that the boundary conditions (11.71) are satisfied. Since 
diam.D m , 1 < m < M, are small, let us write (12.11) as 

M M „ 

u M = u (x) + G(x,x m )Q m + ^2 / [G(x,s) - G(x,x m )]a m (s)ds, (2.2) 

m=l m=l ^ Sm 

where x m G D m is a point inside D m and 

Qm ■= / (Tm(s)ds. (2.3) 



The choice of x m G D m is arbitrary because diamD m < 2a is small. We will 
prove that Q m 7^ 0, give an analytic formula for Q m (formula (12.201) below), and 
approximate the field um in (12.21) by the expression: 



M 

U M = U {x) + ^ G ( X > X m)Qr. 
m=l 



The error of this approximate formula is of order max(|, ka), see estimate (12.71) 
below. Therefore this error tends to zero as a — » since d = 0(a 1//3 ). Let us 
estimate the term 

E m := [G(x,s) -G(x,x m )}a m (s)ds. (2.4) 
By the inequality (I1.19P one gets 

\E m \ < c(Jp + ^j\Q m \, \x-x m \>d^>a. (2.5) 

We will prove below that Q m = 0(a), see formula (I2.20p . and, since \G(x, x m )\ < 
cd" 1 if \x — x m \ > d > 0, one has: 

G(x,x m )Q m \=0(~). (2.6) 



Let us prove that under our assumptions the term E m is much smaller than O(^). 
Using again inequality (11.191) . one gets: 



\E m \ < c(ad~ 2 + kad~ 1 )0{a). 



Therefore, 



,a 2 , a s „,a, 



\E m \<0(- + ka~)<0(~), (2.7) 
because C 1 and a <C d by assumption. So, our claim is verified. Moreover, 

M M 

22 \ E m\ < ^2 \G{x,x m )Q m \ 

m=l m=l 

if \x — x m \ > d 3> a, because M = O(^). 

To find Q m , we use the boundary condition (II. 7p . Let us write u(x) in a 
neighborhood of Sj as 



um(x) = u e (x) + / G(x, ^cr^s)^, |x — Xj| < 2a, (2.8) 
Js 3 

where u e is the effective field acting on the j— th small particle from outside: 

u e (x) := u M {x) - G(x,s)aj(s)ds = ^G(x,x m )Q m + O^). (2.9) 

We neglect the error term O(^) in what follows. From (12.81) and (11.71) one gets: 

= u eN (s) - CM S ) + Ajai 2 ^ - QTjCTj, s G Sj, (2.10) 

where u £ n(s) is the normal derivative of u e at the point s & Sj. One can rewrite 
this equation as: 

Oj = AjUj — 2(jTj<7j — 2(jU e (s) + 2u £ n(s). 
Here the operators Aj and Tj are defined as follows: 

TjOj := j G{s,t)o 3 {t)dt^ j (2.11) 
J Sj J Sj \ \ 

AjCJj := 2 [ a Jt)dt~2 [ ^ -aAt)dt:= Aa h (2.12) 

J J J S j dN s jy } J Sj dN s 47r|s-t| A) 

and we have used the following approximations: 
G(x,y) = g (x,y)[l + 0(\x - y\)], \x - y\ -»• 0; g Q (x,y) :-- 



4tt\x — y\ 



(2.13) 

^^ = ^[l + 0(|x-y| 2 |ln|x- 2/ ||)], |z-y|-0. (2.14) 



Note that (see [9], p. 96, formula (7.21)): 

Aads = - I <7jds (2.15) 



Indeed, 

I ds Ss, wsw^\ G(t)dt = I, dta ' (t) l, ds m- s ^ht\ = - i ^ (()<ii - 

The integral 

/" 9 1 

/ tttt~ ~ — i :ds = —1, t G ib,-, 



is well known in potential theory for surfaces Sj G C 1, . 

Integrating ( I2.10p over Sj, using formula ( 12. 15ft . and the divergence theorem, 
one gets: 

Qj = -CjU e ( x j)\ S j\ ~ Q I ds I -t\ + / AuJy - (2 ' 16) 
The function u e (y) is smooth, so 

/ Au e {y)dy = V j Au e {x j )[l + o{l)}, a -> 0, (2.17) 

where Vj = \Dj\ is the volume of Dj and we have used the smallness of the 
diameter of Dj, that is, the smallness of a. 
Let us write 

(Tj(t)dt f ( \ f ds 



ds / — r = / dtaAt) . 

q. Jq. 4lt\s-t\ Jq. " J q . 4x\S ~ t\ 



<2~ / dt / = ¥4,, Jj := [ I ^*r. (2.18) 

3 Si J q 4tt|s - t| 4^1^-1 3 J Sj Js s \s-t\ 



Here we approximated the continuous on Sj function J s by its mean value 



If Sj is a sphere of radius a, then 

ds 



\s\=a 



= Ana, \t\ = a, (2.19) 
s — t\ 



so in this case equation (12.181) is exact. 

From (12.161) - (I2.18P one finds a formula for Qj\ 



We neglected the term Vj Au e (xj) = 0(a 3 ) which is much smaller than \Q \ \Sj\ = 
0(a) as a — > 0, because \Sj\ = 0(a 2 ) and \Q\ = O(-). The quantity Jj = 0(a 3 ). 
Therefore = 0(1). We choose 

Cj = , (2.21) 

(X 

where H(x) is a continuous function in D, which we can choose as we wish subject 
to the condition ImH < 0, because Im(j < 0. 

If the small particles are all of the same shape and size then \Sj\ = C\a 2 , and 
Jj = c 2 a 3 , where ci, c 2 > are constants independent of j, 1 < j < M, and then 

QJj _H{ Xj )<*. = h{x ^ (222) 



47r|S , j| 4nci 
and 

Cj\Sj\ = H(xj)cia. (2.23) 
Formulas (Q, fl2T20|) . f[2T22|) and fl2T23l imply: 

M - „2„-i 



m=l 



= wo(ac) - ^2 G ( x i Xm "> — TThi — T - Mm ( x ™)' ( 2 - 24 ) 



where \x — x m \ > d ^> a, and we replaced u e (x m ) by Mj\/(x m ) because their 
difference (see (12.91) ) is of order 0(|) <C 1. Indeed 

f c _ a 

\um(x) — u e (x)\ < I \G(x, s)\ \aj(s)\ds < — \Qj\ < c — , \x — Xj\ > d 3> a, 
JSj ' d d 

(2.25) 

where c, c > are some constants independent of a. 
If the assumption (jl.2ip holds, then 



lim G(x, x m ) ^ c f^ UM{Xm)a = f G(x,y) ^-^ u(y)N(y)dy. 
a-o^ l + h(x m ) J D l + h(y) 

(2.26) 

To pass to the limit in (I2.26P one can use the following lemma. 

Lemma 3 Assume that x m G D m , diam.D m < 2a, f is a continuous function in 
D with a possible exception of a point y in a neighborhood of which it is absolutely 
integrable, for example, it admits an estimate \f{y)\ < \ y - yo \b > b < 3, and assume 
that 

lima V 1= / N(x)dx ML) C D (2.27) 

D m CD 
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for any subdomain D C D , where N(x) is a continuous function. Then there 
exists the limit 



M 



f(Vm)a= [ f(y)N(y)dy. (2.28) 

Remark 1 In our case f(y) = G(x,y) u m(v) and f)2.27p is the assump- 

tion dLZU). 

Proof of Lemma [3] Let D = |J p=1 A p , where A p and A 9 do not intersect 
each other, A p is the closure of the domain A p , and linip^oo max p diamA p = 0. 
Choose any point y^i G A p and note that 

sup \f(y {p) )- f(y m )\ <s p ^0 asdiamA p ^0. (2.29) 

£D m ,D m CA p 

Therefore 

M P P 

\imJ2 f(y m )a = lim Vfl £ f(y m ) = X)[/(y«) + 0(e p )} ■ lini La £ 1 

m=l p=l D m cA p p=l D m cA p 

= ^[/(y (p) ) + 0(e p )] / N(y)dy 

p=l Ja p 

p 

P =i 



where linip^oo max p =0. Let P — > oo in (12.301) . Then 

lim Y)[/(y«) + 0(e p )] [JV(y«) + 0(ecj)] |A p | = / f(y)N(y)dy. (2.31) 



In the above argument we assumed that / is continuous in D. If / has an 
integrable singularity at a point x , then we choose a ball B(xq,5 s ) centered 
at x of radius 5 £ such that sup 0<(5<(5e f B r X0 g\ \f(y)\dy < e, where e > is an 
arbitrary small fixed number. Then 



sup / \f(y)\\N(y)\dy<ce, 

0<<5<<5 e Jb(x ,S) 

where c = maXy €D \N(y)\ > is a constant independent of e. Now we apply the 
above argument to the region D\B(x ,5), where / is continuous and get: 

M „ 

lim J2 /W a = / f{y)N{y)dy. (2.32) 

m= l JD\B(x ,S) 
y m £B(xo,8) 

11 



The left side of (12.281) in the case of / having integrable singularity at the point 
Xq and continuous in D\xq is understood as the limit of the expression on the 
left of (E32D as 5 -> 0. This yields (1^28]) . Lemma [3] is proved. ■ 

Passing to the limit M — > oo, or a -> 0, in equation (12.241) and using Lemma 
El one gets 



= Ug (x)-I G(x,y)p(y)u(y)dy } 
Jd 

4tcc 2 %)iV(y) 



p(x) 



c 2 1 + %) 



Applying the operator L = V 2 + k 2 — qo(x) to this equation and using (II. 3p and 
( 11. 4p . one obtains equation ( 11 .23[) . Formulas ( 11.251) and ( 11.261) follow from the 
above equation and from formula (11.121) . 

This concludes the proof of Theorem ■ 

Remark 2 It is possible (and not difficult) to generalize Theorem [2] to the case 
of particles with different shapes. Since this does not lead to an essentially new 
result, we do not go into detail. In [8] one can find analytical formulas for the 
S— matrix for wave scattering by small bodies of arbitrary shapes. 

Proof of Theorem O 

Now we assume ( m = 0, 1 < m < M, which means that all the small particles 
are acoustically hard. In this case equation (12. 10p takes the form 

Oj = AjCTj + 2u eN (s), s G Sj, l<j<M, (2.33) 

where 

u e {x) := uq(x) + / G(x, s)a m (s)ds. (2.34) 

We cannot use approximation (I2.2p because the quantity Q m now is of the same 
order of magnitude as the integral J s [G(x, s)—G(x, x m )} cr m (s)ds, or even smaller 
than this integral. This is established below. While under the assumptions of 
Theorem [2] we had Q m = 0(a), now, under the assumptions of Theorem [31 
we have Q m = 0(ka 2 a 3 ), which is a much smaller quantity than 0(a) because 
fea C 1. To estimate the order of magnitude of Q m , we integrate (I2.33P over Sj 
and use (12.151) . The result is: 

Qj = / u eN ds = I Au e dx ~ Au e (xj)Vj, (2.35) 

JSj JDj 

where Vj is the volume of Dj, and we have used the assumption d ^> a. This 
assumption allows one to claim that u e (x) is practically constant in the domain 
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Dj in the absence of j-th particle. Differentiation with respect to x brings a 
factor k. Since we assume that k > is fixed, this factor is not important for 
our argument, but to make the dimensionality of the term VjAu e clear, we may 
write VjAu e = 0(k 2 a 3 ). This quantity has dimensionality of length since ka is 
dimensionless. 

We now prove that the term E m := J s [G(x,s) — G(x,x m )]a m (s)ds, which 
was neglected under the assumptions of Theorem [21 because it was much smaller 
than \G(x,x m )Q m \, is now, under the assumption ( m = 0, 1 < m < M, of the 
same order of magnitude as \G(x,x m )Q m \, namely 0(/c 2 a 3 <i -1 ), or even larger. 
We have 

[G(x, s) - G(x, x m )]a m (s)ds 
= I V y G(x, a)\ ■ (s — x m ) a m (s)ds, \x — x m \ > d ^> a, (2.36) 

° >~>m 

where we have used the assumption \x — x m \ ^> a and kept the main term in the 
Taylor's expansion of the function G(x, s) — G(x,x m ). 
Recall, that 



s, 



{s-x m ) p a m {s)ds = -V m ^ ) ^^ , (2.37) 

yj Oyj y=x m 



where (3^ is the magnetic polarizability tensor defined in [9], (p. 55, formulas 
(5.13)-(5.15) and p. 62, formula (5.62)), and (s — x m ) p is the p-th component of 
the vector s — x m . Namely, if 



then 



a = Aa — 2Nj, (2.38) 

/ s p a(s)ds = VPtf, (2.39) 
Js 

where V is the volume of the body with boundary S, Nj is the j-th component 
of the exterior unit normal N to S, the role of the point x m from equation (12.371) 
is played by the origin, which is located inside S, and the role of S m is played by 
S. Equation (I2.33j) with j = m can be written as 

a m = A m a m - 2N 3 (-^^) • (2.40) 



Compare (E3QD and (12381) and get (l2~37l) . 



Formulas for the tensor j3 P j = oc P j(j) for bodies of arbitrary shapes are 

7=-l 

derived in [9], p. 55, formula (5.15), so one may consider the tensor (5 P j known for 
bodies of arbitrary shapes. The parameter 7 = £i ~ €e , where is the dielectric 



13 



permittivity of the body and e e is the dielectric permittivity of the surrounding 
medium. The case 7 = — 1 occurs when q = 0. This is the case, for example, 
in the problem of calculation the magnetic dipole moment of a superconductor 
placed in a homogeneous magnetic field: in the superconductor the magnetic 
induction vector B = 0, which means that the magnetic permeability ^ of such 
body is zero, = 0, see [2]. That is why the tensor [3 P j is called magnetic 
polarizability tensor in [9]. 

From (1X361) and f[2lITj) it follows that 



[G(x,s) - G(x,x m )]a m (s)ds ' ({ ' J> 



y=x„ 



V (m) 



y=x„ 

(2.41) 

where one sums up over the repeated indices p, j, but nor over m. The quantity on 
the right of (12.411) is of the order 0(k 2 a 3 d~ l ) if kd > 1, that is, of the same order 
as \G(x, x m )Q m \, provided that \x — x m \ > d ^> a, and it is of the order 0{ka 3 d~ 2 ) 
if kd < 1. Indeed, fyf = 0(1), V m = 0(a 3 ), and \V v G(x,y)\ < cmax (J, 
Let us prove the estimate 

k 1 

\V y G(x, y) J < cmax ( — , — J for \x — y\ > d ^> a, 

where c > is a constant independent of d. 
We have 

G(x,y) = g(x,y) - / g(x, z)q Q (z)G(z, y)dz := g - TG, 
Jd 

where T is compact as an operator in L p (D),p > 1 under our assumptions, 
namely, D C R 3 is a bounded domain, qo(x) is a bounded piecewise-continuous 
function. From this equation we get 

V y G{x,y)=V y g{x,y)-TV y G. 

Clearly, 

V7 ft 1 ^ y~ x 
\/ y g = g(ik 



so 



\x — y\ \x — y\ 



\V y g(x,y)\ < 2max (A, J^) = i-max (£, 1), |* - y| > d > 0. 



Thus 



\V y G{x,y)\ < \V y g(x,y)\[l+c / ——\V y G(z,y)\dz .., ^ ~ *!' 77 ] := |V^|(l+cJ) 
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where 

I:= / \V y G(z,x)\ - ? =f=^ <c 
>D \ x — z \ a/1 + K l \x — yr jd 



dz 


x — 


y\ 


2 


\x — z\ 


V 1 + k2 \ 


X 


- y\ 2 



dz 


x-y\ 


2 


\z - y\ 


2 


x — z\ 


y/1 + k 2 d 2 



One has 

dz 



D \z — y\ 2 \x — z\ 



< c In \x — y\ 



and 

sup | In \x — y\ | \x — y\ 2 < c, 

x,y&D 

where c = c(D) is a constant. Therefore 

c 

I < - < c, 

~ VI + k 2 d? ~ 

and II 1 

\V v G(x, y)\ < cmax (—, —) . < cmax (-, — ), 

1 v v ' y;| - v d' d 2 ' y/i + k 2 d 2 ~ y d d 2> 

as claimed. 

H | > -h, i-e. fed > 1, then \V y G(x,y)\ <c\, \x - y\ > d > 0. 

If | < ^, i.e. kd < 1, then \V y G(x, y) \ < Jr, |x — y| > d > 0. 

Therefore, the right side of (12.411) is O(^-j-) if fc<i > 1, in which case it is of 
the same order as the term G(x, x m )Q m . If kd < 1, then the right side of (I2.4ip is 
O(^pr), m which case it may become larger than the term G(x, x m )Q m because 
the ratio — ^ > 1 provided that kd < 1. 

Writing the field (12.1 1) in the form 

M M , 

u M (x) = u {x) + ^2G(x,x m )Q m + ^2 / [G(x,s)-G(x,x m )]a m (s)ds (2.42) 
and using formulas (12.351) and (12.411) . one gets: 

u M {x) = u (x) + ^ G ( x > x m )Au e (x m ) V m -^2 dG ^ ,Xm ^ dUe ( Xm ^ V m p$\x m ), 

m=l m=l Vp Vj 

(2.43) 

and over the repeated indices p,j one sums up. 

Let a — > 0, M — > oo. We want to give sufficient conditions for passing to this 
limit in fl2^3l) . 

Lemma 4 ^ss^me i/iat for any subdomain D <Z D the following limits exist: 

lim V V m (3^\x m ) = [ P PJ (y)v(y)dy, (2.44) 

lim Y J V m =( u{y)dy. (2.45) 

a^O y 5 

OrnCf 
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D m CD 



Then the limiting form of equation (I2.43P is: 

U{x) = u {x)+ [ G{x,y)AU{y)v{y)dy 

J D 

dG(x,y) dU(y) 



Id dy P dyj 
where one sums up over the repeated indices p, j . 



f3 pj (y)u(y)dy, (2.46) 



Remark 3 If one assumes that u(y) vanishes near the boundary S of D and 
integrates the last integral in ( 12.461) by parts, one gets 

U(x)=u (x) + / G(x,y){AU(y)u(y) + V ^(^Mp pj (y) u(y))\dy. 

Jd 1 p j=l °Vp v °Vj J } 

(2.47) 

Applying the operator L to both sides of (12.471) and using (11.41) one gets: 

uu + u(y)Au(x) +J2w~ {nr^Mv) "(f)) = °- ( 2 - 47 ') 

9 j =1 °Dp v °Dj 



Remark 4 If all the small particles are identical, then V m = c%a , where the 
positive constant C3 does not depend on m, and (3^ = (3 P j. Then 



lim V V m = lim Cl a 3 V 1 = \im[ Cl a 3 N(D)\, (2.48) 



D m CD D m CD 



where N(D) is the number of small particles in the domain D. For the limit 
( I2.48P to exist it is sufficient that 

N(f)) = f D y(y)dy 

Cia A 

where u{y) > is a continuous function, and the limit in (I2.48[) is equal to 

I D v(y)dy. 

One can write (12.491) as 

N(y)dy=^ldy. (2.50) 
cicr 

In contrast to Theorem 2, where M = O(^), we now have M = O(^). 
Similarly, 

lim V V m ^f(x,) = \im[ Cl a 3 (3 PJ N(D)}=P PJ [ u(y)dy. (2.51) 

^— ' J a^O J rj 

DmCD 
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We gave in this Remark some practically realizable sufficient conditions for the 
existence of the limits (12.440 and (12.450 . 

Let us verify that if the limits (12.441) - (12.450 exist, then the limit of the right 
side of equation (12.430 exists, and, denoting this limit by 

U.[x) = lim-UM^), 
one obtains the limiting form of equation (12.430 : 



U{x) = u (x)+ / G(x,y)AU(y)is(y)dy 

J D 

f dG(x,y) dU{y) ( w . 

-qT Qj-MVMV)*, (2-52) 



which is equation (12.460 . 

We took into account that 



lim u e (x) = U(x) 



This is so because, as a — > 0, the input of a single particle into the field U{x) 
tends to zero. 

Let us verify the existence of the limit of the right side of equation (12.430 . We 
use, as in the proof of Lemma El a representation of D of the form D = |J x A p , 
and assume that 

lim max diamAp = 0. (2.53) 

P-+oo l<p<P 

Then 

M 



lim V] G(x, y m )Au e {x m )V n 

a— >0 ' ' 



m=l 





p 

= V"lim y~] G(x,x m )Au e (x n 

P=l D m CAp 






P 

= J2G(x,y^)Au e (y^)(l + e p ) 
P =i 


lim V m 

-Dm CA P 




p 

= Y,G(x,y {p) )Au e (y^)(l + e p ) 
P =i 


v 2 (y (p) )(l + e' p )\A p \ 


where 


lim max( \e„\ + \e'\) 

P^oo p 


= 0. 



(2.54) 



(2.55) 
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Let P -> oo in f[2~54l) and use (123511 to get 

p 



lim £ G(x, yW) A^yW) | A p | (1 + £ p + ej + e^) 

p=l 

= [ G(x,y)AU(y)u(y)dy. (2.56) 
Jd 

We have replaced u e (y^) in the limit P — > oo by W(y), because 

U(y)-u e {y) = G(s,t)a(t)dt = o(l) as a -> 0. (2.57) 



From (12.56P and (12. 54j) one gets: 

a/ „ 
lim VG(x,x m ) A Me (x m )\/ m = / G(x,y)AU(y)is(y)dy. (2.58) 

a_>0 — /r> 
m=l 

The singular points x = y G D of G(x, y) are treated as in the proof of Theorem 

m 

The function \G(x,y)\ < c\x — y\~ l as |x — y\ — >• 0, so |G(x,y)| e L X {D) as a 
function of y for any x £ D. 

Similar arguments, applied to the last sum in ( 12.4=31) . lead to the formula 



^ dG(x,x m ) du e (x m ) (m) r dG(x } y) dU , )() , 

-OA d y p d y. P"P3 ^m) J d 

(2.59) 

where one sums up over the repeated indices p,j. 

Theorem [3] is proved. ■ 

In Section H] we discuss the compatibility of the condition d ^> a and the 
existence of the limit ( 11 .281) . 

3 Auxiliary results 

In this Section we prove Theorem [1] and Lemmas [H, EJ 

Proof of Lemma [1] Let us start with the following observations: 



e 



t — y| — \x — y\ \ < \t — y — (x — y) \ = \t — x\ < a, (3.1) 

sup |e is -l|<a, (3.2) 

a<s<a 

ik\t- y \ _ e ik\x- y u = | e ifc(|t-i/|-|*-i/l) -\\<ka, (3.3) 
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where the last inequality follows from (13.21) . 
One has 

■ , I \x - y\e ik \ l -y\ - \t - y\e ik \ x ~ y \\ \\x — y\ — \t - y\\ 

M'-d-^l = \*-v\\t-v\ \*-y\\t-y\ + 

|£ _ y|| e ifc l*-y| _ e ik\x-y\ 

+ \x-y\\t-y\ (3 - 4) 
a ka a ka , a , ka 

~ \x-y\\t-y\ + \x-y\ ~ <P(l-±) + T ~ ^d?' + T' 

Lemma [1] is proved. ■ 

Proof of Lemma [2] Let us start with the equation: 



G(?,y)=g(x,y)- g(x, z)q (z)G(z,y)dz, (3.5) 
Jd 

where go is defined in (11.21) . From (I3.5P one gets: 

\G(t,y)-G(x,y)\ < \g(t, y) - g(x, y)\ + I / [g(t, z) - g(x, z)]q (z)G(z, y)dz\ 



D 



a JvCl f dz 

s °^ )+ ^ +c j D w - z) - 3[x - z) ^[ < 3 - 6) 



Here we have used Lemma [T] and the estimates 



sup \q (z)\ < c 4 , \G(z,y)\<c 5 \z-y\ 1 , (3.7) 
zeD 



where C4, C5 > are some constants. 
Let us estimate the integral 

I ■= [ \g(t,z) - g(x,z)\ 



D 



z-y\ 



\g(t,z) - g(x,z)\dz f \g(t,z) - g(x,z)\dz 



x-z\>f,zeD \ Z — U\ J\x~z\>±,zeD \ z — U\ 

:= h + h- (3.8) 
By Lemma [U which is applied to I\ with d replaced by |, one gets 

_ ^ / a ka, f dz , a ka, ,„ _ 

''^k + T^U^^k + T)- (3 ' 9) 

Here and below we do not write z G D under the integration sign to simplify the 
notations. 
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Let us estimate I2: 



x-z\<f 



dz\ 






x — z\ 




\t-z\\ 




z - 


- 2/11* - 


- z\ 


\x — z\ 





4:71 

One has 

\e tklt - zl \x - z\ - \t - z\e iklx ~ zl \ < \\x-z\ - \t - z\ \ + \t - z\ |e lfc| *~ z| - e iklx ~ zl \ 

< \x-t\ + \t-z\k\\t-z\ - \x-z\\ (3.11) 

< \x-t\ + k\t-z\\t-x\. 

Thus, with \x — y\ > d ^> a and \t — x\ < a, one has: 
1 



\x-z\< 



dz(\x -t\ + k\t - 


-z\\t- 


- x\ 


) 




z - 2/11* - A 


\x — z\ 





\t — x\ ( f dz , f dz 



4?r ^ J\*-z\<$ \ z - v\ \t - A \ x - A J\x-z\<i \z-y\\x-z 

< ca{± + \). (3.12) 

From f)3.9p and (13.121) the estimate (j!.19p follows. Lemma [2] is proved. ■ 

Proof of Theorem [TJ Let us first prove that if conditions (11.141) hold, then 
problem (11.51) - (11.71) has at most one solution. It is sufficient to prove that the 
homogeneous problem 

M 

(V 2 + k 2 - q )u = in M 3 \ [j D m , (3.13) 

m=l 

du 1 1 

— iku-o(-), u-O(-), r = \x\ — > 00, (3.14) 

un = CmU on S m , 1 < m < M, (3.15) 

has only the trivial solution if conditions ( 11.141) hold. 
Taking complex conjugate of ( 13.131) - (13. 15ft one gets: 

M 

(V 2 + k 2 - %{x))u = in M 3 \ [J D m , (3.16) 

m=l 

du I 1 

— — \-iku — o(-), u-O(-), r = \x\ — ► 00, (3.17) 
un = C m u on S m , 1 < m < M. (3.18) 

Multiply ( I3.13P by u, (I3.16P by u, subtract from the first equation the second 
one, and integrate over the region (IR 3 \ (Jm=i D m ) H -B/? := -D_r, where is the 
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ball centered at the origin of radius R. Using Green's formula, one gets: 

= / [uV 2 u — uV 2 u — (qo — %)\u\ 2 ]dx 
Jd r 

f f du du 

= —2i / lmqo(x)\u\ dx + / (u— u—)ds 

JDu J\x\=r or or 



E' f , du du x , . 
, ^m- u m^- (319) 

m=l J ° m 



Using (13.171) and (I3.18P one rewrites (13 . 1 9[) as follows: 
= — 2% / lmq Q (x)\u\ 2 dx + 2ik / \u\ 2 ds + o(l) — 2i^^ / Im^-u^ds. 

JDr J\x\=R m=1 JSm 

(3.20) 

Letting R — ► oo, taking into account that qo(x) = in D' = M. 3 \D, and one gets: 



M 



< f Lna.WH'dx + V / Im Cm |«|^ - tlimsup / |«|»*. 

(3.21) 

Since all the terms on the right side of this relation are non-positive by the 
assumptions ( 11.141) . it follows that 

limsup / |w| 2 <is = 0. 

R-^co J\x\=R 

This implies that u = (see, [8], p. 231). 

Thus, uniqueness of the solution to problem (II. 5p - 1.7 is proved. 

Let us prove the existence of the solution to (I1.5P - (11.71) of the form (ll.lOp . 
The existence of the solution of the form (ll.lOp will be established if one proves 
the existence of a m , 1 < m < M, such that boundary condition (ll.7p is satisfied: 

U eN ~ CjUe + ~ ~ QTjCTj =0, 1 < j < M. (3.22) 

Here u e , which depends on j, is defined by the formula: 

/ G(x, s)<Tj(s)ds = u + / G(x, s)a m (s)ds. (3.23) 



u P := u — 



Yl,A 



Under our assumptions S m G C > uniformly with respect to m. Therefore equa- 
tion (13.221) is of Fredholm type in the space £ 2 (IJm=i ^m)- The corresponding 
homogeneous equation, i.e., the equation with m = 0, cannot have a nontrivial 
solution because such a solution would generate by formula (ll.lOp with % = 0a 
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function um(x) = Y^ m =i Is ^( x > ^^mi^ds, which would solve the homogeneous 
problem (jl.5p - (jl.7p . We have already proved that such a function has to be zero 
in R 3 \ IJm=i ^m- Thus, WmI = 0, 1 < m < M, and Um solves the problem: 

Lqu m = in D m , um\ = 0, 1 < m < M. (3.24) 

If diam.D, m < 2a is sufficiently small, then problem (13.241) has only the trivial 
solution for every m, 1 < m < M. Therefore Um = in D m and in R 3 \ (J^ =1 D m . 
Therefore, by the formula for the jumps of the normal derivatives of the single 
layer potential, 

dN \s m dN \s m 

we conclude that a m = 0, 1 < m < M . This implies the existence of the solution 
to problem flTJS} - (HHD of the form (TTTT0|) . 

Theorem [1] is proved. ■ 



Let us return to the assumptions of Theorem 2, namely, 

( m = O(-), aN(A b (y)) = N(y)\A b (y)\(l + o(l)), 
a 

where A b (y) is a cube, centered at the point y with the side b > 0, and o(l) is 
related to the limiting process 6 — > 0. 

Under these assumptions let us establish an estimate for the function vm '■= 
um — Uq, which is uniform with respect to M — > oo, or a —>■ 0. From this estimate 
it follows that Vm converges, as a — > 0, in L 2 (M 3 , (1 + |a;|) _1_7 ), where 7 > is 
an arbitrary fixed constant. The function vm satisfies the radiation condition at 
infinity. The function m G Hf oc {R 3 ) solves the equation L uo = in R 3 . 

Let D e := R 3 \ U^ =1 D m and S' := u£f =1 S TO . Let 

\\ v \\ := ([ \v(x)\ 2 {l + \x\y l -Ux^j 112 , \\\v\\\ = J2([ {\v N \ 2 + \v\ 2 )d s y /2 . 

The estimate we wish to prove is: 

\\ v m\\ < c|||m |||. (3.25) 

Here and below c > stand for various constants independent of a. 
Let us outline the proof of inequality (I3.25p . 

Step 1. If M = O(J), then the right side of IEE2511 is bounded as a — > 0. 

Indeed, the number of small particles is M = 0(-) and uq is Hf oc {R 3 ), so that 
m and uon are bounded in L 2 (S m ) uniformly with respect to m, 1 < m < M. 
Thus, 

C f f \ ^/^ c 

||ko|||<- max ( / (|«oiv| 2 + |mo| 2 )^ ) < -\S m \ 1/2 < c, 
a Km<M \ q ) 



a 
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where c > stand for various constants independent of a. 

Step 2. If the inequality (I3.25j) is false, then there is a sequence u^, \ \ \u^\\ \ = 
1, such that \\v^ || := ||f^|| > n. 

Define w (n) := p^y. Then 

||w (n) || = l. (3.26) 

From the weak compactness of bounded sets in L 2 , it follows, that one may select 
a subsequence, denoted again w^ n \ such that converges weakly in Lf oc (D') 
to a function w. The function w^ n ' solves the problem: 

L w {n) = in D e , 
-( m w^ = (( m u { ^ -u^)/\\v^\\ on S m ,l<m<M } (3.27) 

and satisfies the radiation condition. 

It follows from fl3T27D that ||VV n )|| < c, so \\w^\\ H 2 oc{De) < c, where H 2 oc (D e ) 
is the Sobolev space. Thus, one may assume, using the compactness of the em- 
bedding from Hf oc into L 2 oc , that converges to w strongly in Lf oc (D e ). This 
and equation (I3.27P imply that converges to w strongly in H 2 oc (D e ), so that 
w solves equation (I3.27p . satisfies the radiation condition and the homogeneous 
boundary condition (13.271) . that is, wn — ( m w = on S m . Therefore, by al- 
ready proved uniqueness theorem (see the proof of Theorem 1), we conclude that 
w = 0. The terms Uq /\ \v^\ | and u^/\ \v^\ | tend to zero as n — > oo, be- 
cause ||f (ri * ) || > n. Therefore, the limiting function w satisfies the homogeneous 
boundary condition wn = CmW on S m , 1 < m < M. 

Let us prove that |u/ n )(:r)| < -A, \x\ > R, where R > is sufficiently large 
and c > does not depend on n. 

For one has a representation by the Green formula in the region \x\ > R, 
where R > is large enough, so that the ball Br := {x : \x\ < R} contains D. 
Namely 

w V>(x)= I (wl n) g(x,s)-g r (x,s)w {n) )ds, \x\ > R, (3.28) 

J\s\=R 

where the derivatives with respect to r are the derivatives along the normal to 
the sphere Sr := {s : \s\ = R}, and g is defined in (11.151) . It follows from (13.281) 
that |u7^(x)| < ■At for |x| > R, where c > is a constant independent of n, 

because local converegence in H 2 implies that the L^S^-norms of w^ n ' and of 
Wr are bounded uniformly with respect to n. 
Therefore 

lim \\w {n) - u; 1 1 = 0, (3.29) 
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because on compact sets lim n ^oo \ \w^ n ' — w\\ H 2 ( D n = 0, and near infinity the 
inequality Itw^^x)! < ^ implies that 

/ |w (n) (x)| 2 (l + Ixl)" 1 " 7 = 0(iT 7 ) -»■ 0, i? -> oo, 

J{x:\x\>R} 

so that (I3.29f) holds. Because of the uniqueness of the limit, not only a subse- 
quence of but the sequence itself converges to w as n — > oo. 

This leads to a contradiction, because w = and f!3.26|) together with (13.291) 
imply | \w\\ = 1. 

This contradiction proves inequality (13.251) . 

From inequality (13. 25ft and Step 1 one concludes that that um contains a 
weakly convergent in Lf oc (D e ) subsequence. By the arguments, similar to the 
given above, this subsequence converges in L 2 (M 3 , (1 + |x|) -1-7 ). Its limit solves 
equation (11.231) . 

The relation M = O(-) plays an important role in our proof of Theorem 2 
and in Step 1 in the above argument. 



4 Application to creating smart materials 

Let us ask the following question: can one make a material with a desired re- 
fraction coefficient n(x) in a bounded domain D C M 3 , filled by a material with 
a known refraction coefficient n (x), for example n (x) = n = const in D, by 
embedding into D a number of small particles, each of which is defined by its 
shape and its boundary impedance? 

Consider first the particles satisfying the assumptions of Theorem [2j More 
specifically, suppose that all the particles are balls of the same radius a. In this 
case ^ 

\S m \ = Ana 2 , \J m \ = / = 167r 2 a 3 , 

J \s\=a J \t\=a 

SO 

n 47TC? 



\s-t\ 



c\ = An, c 2 = 16n , = An, 

c 2 



and formula (jl.24p yields 



1 + h(x) 



where h(x) is defined by the choice of the boundary impedances by formula (ll.20p : 

CW = ^, (4.2) 
a 

and N(x) is defined by formula (ll.2ip . 
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If the original refraction coefficient is uq(x), then the corresponding potential 
is qo(x) = k 2 [I — n (x)} by formula (II. 2p . If the desired refraction coefficient in 
D is n(x), then the corresponding potential is q(x) = k 2 [l — n(x)], so 

p(x) = q(x) — qo(x) = k 2 [no(x) — n(x)]. (4.3) 

To create a material with the desired refraction coefficient n{x) it is sufficient 
to choose N(x) and h(x) so that (14.11) holds with p(x) defined in (14.31) . If the 
new material with the refraction coefficient n(x) has some absorption, that is, 
lmn(x) > 0, and lmn = 0, then lmp(x) < 0. Let us prove that any function 
p(x) in D with Imp < 0, can be obtained (in many ways, non-uniquely) by 
formula (14.11) with some choices of a nonnegative function N(x) and a function 
h(x) with Im h < 0. 

Let p(x) = pi(x) + ip 2 (x), p 2 (x) < 0, and h(x) = hi(x) + ih 2 (x), h 2 (x) < 0. 
Assume that p(x) is given. Then (14.11) implies 

(hi + iho) (1 + hi — ihn) AT/ N , s 

Thus 

» ^ 4 * (TttStI ^ 4 'iiwn^' (4 - 5) 

There are many choices of the three functions: N(x) > 0, h 2 (x) < and a real- 
valued function hi(x) such that relations (14. 5p hold. For example, if p± > and 
£>2 7^ 0, then one can choose 

ftlW= 0, JV W = ?M±fiM. (4.6) 

P2(a;) 47rpi(x) 

It is a simple matter to check that relations (14.51) hold with the choice (14. 6p . 
Since one has three functions hi(x), h 2 (x) < and N(x) > to satisfy two 
equations (14.51) with p 2 (x) < 0, there are many ways to do this. A particular 
choice of h(x) = h\(x) + ih 2 (x) and N(x) > yields the surface impedance 
((x) of the particles to be embedded around each point x G D, ((x) = by 
formula (14. 2p . and the number of particles per unit volume around the point x, 
namely, by formula (jl.2ip this number is so that the number of particles to 
be embedded in the volume dx around point x is equal to dx. The smallest 
distance d between the embedded particles should satisfy the inequality d a. 
One may try to take practically d > 10a. 

Example 1 Suppose that the elementary subdomain A p , used in the proof of 
Lemma [31 is a cube with the side b 3> d, x G A p . Let, for example b = 10~ 2 cm, 
d = 10 _3 cm, a = 10 _5 cm. Then there are (^) 3 = 10 3 small particles in A p 
around a point x, the center of A p . The function N(x) in A p in this example 
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is found from the formula = 10 3 (use (11.211) with D = A p ), so N(x) = 

10 -5 • 10 3 • 10 6 = 10 4 . The number of small particles, embedded in the cube 
A p around point x, the center of this cube, is 10 3 in this example. The relative 
volume of these particles in A p is 10 3 • 1 7rl0~ 15 • 10 6 = 4.18 • 10~ 6 , so it is quite 
small, which is in full agreement with our theory. 

The assumption (II. 8p . specifically, d 3> a, is compatible with the requirement 
f lL2T]) . Indeed, if one denotes by N(D) the left side of OB . then N(D) = 0(±) 
for any D C D. 

Let us assume that D is a unit cube, and denote by N(D) the left side of 
(ll.2ip . The assumption d 3> a implies that the number N(D) of particles in 
D is O(Jr). These relat ions are compatible if and only if OQ) = O(Jj), i.e., 
d = 0(a 1 ^). Therefore, it is possible to have a — > 0, | — > and equation (11.211) 
satisfied. 

Let us discuss the new material properties, specifically, anisotropy, when 
acoustically hard particles are embedded in the domain D, and the assumptions 
of Theorem [3] are valid. The physical situation is now quite different from the 
one in Theorem [2J From the physical point of view one can anticipate the drastic 
difference because the wave scattering by one small acoustically soft particle of 
the characteristic size a is isotropic and the scattering amplitude is of order a, 
while the wave scattering by a small acoustically hard particle is anisotropic and 
the corresponding scattering amplitude is of order k 2 a 3 , (see [9], chapter 7). We 
assume that ka <C 1, say ka < 0.1, so that the quantity k 2 a 3 = (ka) 2 a is 100 
times less than a. 

Example 2 Let us assume again that the small particles are all balls of the same 
radius a. Then 

V m = -ira 3 , is(y)\A p \ = -na 3 N(A p ), 
where iV(A p ) is the number of small particles in a small cube A p centered at the 

3 

point y. If b is the size of the edge of the cube A p , then v(y) = 4.18 p- N(A P ), 
where 4.18 is an approximate value of -y . The magnetic polarizability tensor f3 P j 
of a ball of radius a is j3 P j — — | 5 P j , while the electric polarizability tensor of a 
perfectly conducting ball is 35^, where 

5 P j = 

These values differ by the factor 4n from the values in |2] because we use the 
formula if = for the potential of a point charge, while in [5] this potential is 

A . In our example [3 P j does not depend on m. Therefore the limit (11.271) exists 
if the limit (11.281) exists. The limit (I1.28P exists if and only if the following limit 

exists: 

-^lim a 3 V 1= / v(y)dy, (4.7) 

D m CD 
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where u(y) is the function defined in (11.281) . Thus, in contrast to Example (TJ 
where N(A p ) = now we have N(A P ) =0(£). The relat ive volume of the 

small particles in Example [2] is not negligible and does not go to zero as a — ► 0, 
in contrast to Example [U 

Let us discuss the compatibility of the condition d 3> a and the existence 
of the limits (I1.27P and (11.281) . If the condition d ^> a is compatible with the 
existence of the limit (11.281) . then it is compatible with the existence of the limit 
(fL27j) . If the limit (OHD exists, then a 3 N(D) = 0(1), so N(D) = 0(a~ 3 ). On 
the other hand, N(D) = 0(d~ 3 ). These relations, in general, are not compatible 
because d 3> a. Let us argue more precisely. Let D = A p , where A p is a cube 
with the edge of size b. Let us assume that the small particles in A p are identical 
and their characteristic size is a. If (I1.28P holds, where v(y) is continuous, and if b 
is small, then the right side of (I1.28P equals to z/(y)6 3 , y <E A p = D. The left side 
of (TL281) equals to c 3 a 3 N{A P ). Thus iV(A p ) = i v(y) § . On the other hand, 

N(A p ) = ^ , prvided that one assumes that the centers of the small particles 
are at the uniform grid, so that there are | centers on the segment of length b. 

If ^ v(y) ^ = ^ , then | = (^^■) 1 ^ 3 - Therefore the condition d 3> a is satisfied 

only if (^-) ^ 1, say (^-)^ 3 < 0.1. The number c 3 depends on the shape 
of the particle. If the particles are balls of radius a, then c 3 = 4.18. Therefore 
u(y) < 4.1(T 3 . 

The conclusion is: 

The condition d ^> a is compatible with the existence of the limit ( I1.28P only 
if the function u(y) in ( I1.28P is sufficiently small. 

In general, equation (ll.30p cannot be reduced to a local differential equation 
for U{x). However, if v{y) is small, one may use perturbation theory to study 
equation (11.301) . However, under an additional assumption, reasonable from the 
physical point of view, one can reduce integral-differential equation (ll.30p to a 
differential equation. Namely, let us assume that v(y) is a continuously differen- 
tiable function in D which vanishes near the boundary S. 

Under this assumption one can integrate by parts the last integral in (11.301) 
and get: 

U(x)=u (x)+ [ G(x,y)\AU(y)p(y) + V A /3 pj (y)v(yj)] . (4.8) 

Jd l p j=1 cy P v oyj / J 

Let us apply the operator L Q = V 2 + k 2 — qo(x) to (14.81) and use (11.41) to get: 

[V 2 + k 2 - q (y)]U + v(y)V 2 U(x) + £ A (^A My) u{y ^ = , (4.9) 

P ,j=i Vp y 
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where U[x) satisfies the radiation condition of the type (11.61) . This is an elliptic 
equation and the perturbation V of the operator L is: 



This perturbation is the sum of the terms with positive small coefficient v(y) in 
front of the second derivatives of U and a term with the first order derivatives of 
U: 



If both v{x) and Vv(x) are small, this equation can be studied by perturbation 
methods. The physical effect on the properties of the new material, created 
by embedding into D small acoustically hard particles, consists in appearing of 
anisotropy in the new material: the propagation of waves is described by the 
integral-differential equation (14. 8 p or (under the additional assumption on v(y), 
namely: u(y) vanishes near the boundary S of D) by the differential equation 
(14. 9 p with variable coefficients in front of the senior (second order) derivatives 
and the terms with the first order derivatives. 

The role of the compatibility of the assumption d 3> a and of the assumption 
(I1.28P is quite important. Although passing to the limit a — > 0, justified in the 
proof of Theorem [3], is based on the assumptions (11.271) and (11.281) . but without 
the assumption d^$> a one cannot expect, in general, that the effective field u e (x), 
acting on any single particle, is practically constant on the distances of the order 
2a. This physical assumption is important for our theory. 

From the mathematical point of view, if u(x) is not sufficiently small, then the 
existence of the unique solution to equation (14.81) or of the solution to equation 
(14. 9p . satisfying the radiation condition, is not guaranteed. 

If, on the other hand, the quantity 



that is, this quantity is sufficiently small, then one can argue that the norm of 
the integral operator in (14. 8p in L 2 (D) is small, so that equation (14.81) has a 
unique solution in L 2 (D). This solution admits a natural extension to the whole 
space M 3 by the right side of (14.81) because u(y) vanishes outside D. Since G(x, y) 
satisfies the radiation condition, the solution to (14. 8 p also satisfies this condition. 
Without the assumption that |^(a;)| + |Vz/(x)| is sufficiently small, one cannot 
use the above argument. 

With this assumption one may solve equation (14.81) by iterations and find in 
this way an approximate solution to this equation. The first iteration yields the 




(4.10) 




(4.11) 



sup (\v(x)\ + \Vv(x)\) < 1, 
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following approximate solution to equation fl4.8[) : 



U{x) = u (x) + / G(x,y) Au (y)iy(y) + ^2 



d / dupjy) 
dy p \ dyj 



Ppj(y)v(y))\dy- 

(4.12) 



Formula (14.121) gives the correction to the solution Uq(x) of the unperturbed 
scattering problem, i.e., the scattering problem in the absence of small bodies. 
Since one has 



In [3], Chapter 3, Section 3, the Neumann problem for the Helmholtz equation 
with n (x) = 1 was studied in the domain, similar to the one in equation (11.51) 
and it was proved under the assumptions used in [3], that the main term of the 
asymptotics of the solution, as the relative volume of the particles tends to zero, 
is the incident field, while the next term is proportional to this relative volume. 
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